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Abstract: Relativistic prescription is used to study the slow rotation of stars composed by self- 
gravitating bosons and fermions (fermions may be considered as neutrons). Previous results demand 
that purely boson stars are unable to display slow rotation, if one uses relativistic prescription with 
classical scalar fields. In contrast to this, the present work shows that a combined boson-fermion 
star in its ground-state can rotate. Their structure and stability are analysed under slow rotation 
approximations. 
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Stars composed by neutrons (fermions) has been first suggested by Oppenheimer and Volkoff pL After that, many 
theoretical studies have arisen on this theme and some years later many stellar objects has been identified with neutron 
stars. Inside neutron stars, fermions are attracted each other gravitationally but the configuration is prevented from 
collapsing via Pauli exclusion principle pressure. 

Another possibility for stars configurations has been suggested Q independently by Kaup, and Ruffini and Bonaz- 
1 zola. Boson stars are compact configurations made with bosons gravitationally bounded. Bosons are described via 
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a scalar field, and large interest has been shown on this subject since properties of that kind of configurations can 
be obtained directly from the Lagrangian treatment, without the need to define an equation of state, as it is usually 
£SJ I done to ordinary stars (see also Breit et al ||). Boson stars are prevented from collapsing gravitationally due to a 
. pressure that arises from the Heisenberg uncertainty principle. Boson stars formation by Jeans instability has been 
studied, e.g., by Grasso (Tt]] and its stability has been addressed by Gleiser and by Torres ||. Some reviews |J] 
present the relativistic and astrophysical interest on this subject, properties and detection expectations. 

Since in the primordial gas possibly bosons and fermions coexisted, one could expect stable configurations of self- 
gravitating bosons and fermions formed. A model boson fermion star has been first proposed by Henriques, Liddle 
and Moorhouse || . Some studies on stability || and on configurations with interactions between bosons and fermions 
|l5| has been performed. 

5h ■ Boson star rotation has been first addressed by Kobayashi et al ]Tl| ] . They used a perturbative ap pro ach and have 
jpf^' shown that boson stars could not display slow rotation. But, as shown by Silveira and de Sousa |l2| , considering 
. the scalar field quantum nature with axial symmetry it is possible to obtain stable configurations with I ^ 0, using 
Newtonian approach. After that, many works have contributed to the comprehension of the properties of spinning 
boson stars, as Ryan |D| which shows the large self-interaction case, and Yoshida and Eriguchi Jl4j studied static 
■ relativistic axisymmetric solutions. 

But it seems there is no explicit proof that boson fermion stars can display rotation. Hence, in this article we 
present our results for relativistic slowly rotating stars composed with bosons and fermions. In section II we present 
the construction of a boson star in II. A and of a boson fermion star in section II. B. In section III we present the 
construction of a rotating fermion star in III. A, and an extension for studying rotating boson fermion stars in III.B. 
Results are shown in section IV. 
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II. NON-ROTATING BOSON-FERMION STARS 

The study of composed boson-fermion stars has been first addressed by Henriques et al || . In this section we present 
a summary for boson stars model and its enhance for boson-fermion stars by introducing fermions with Chandrasekhar 
perfect fluid model. 



* Electronic-addresses: claudio@unb.br, claudio@iccmp.br, desousa@ucb.br 
t Electronic- address : vandalrs@tcu .gov .br 



1 



A. Boson Stars 



A non-rotating configuration is spherically symmetric and one can build the metric into the form: 

ds 2 = -B(r)dr 2 + A{r)dr 2 + r 2 d6 2 + r 2 sin 2 6dip 2 (1) 

Our sign and conventions are the same as those used by Liddle and Madsen [Q. The Lagrangian density for the 
complex scalar field with mass m, is: 

R 



C 16nG 

and, as usual for boson stars, we use the ansatz: 



$(r, t) = <p(r)e 



(2) 



(3) 



The complex field <j>(r) may be written as a sum of real and imaginary parts. In absence of interaction between 
bosons and fermions, real and imaginary parts generate similar sets of equations. This way one can consider <f>(r) 
purely real; this is not possible if fermions and bosons interact Jl5| since <f>i (r) and (f>2 (r) present different behaviour 
(considering <p(r) — <pi(r) + ifafr) above). 

For simplicity we are not considering scalar field self-interaction term, A$ 4 , and we are not considering boson- 
fermion interaction although it should give more realistic results fTo)] . This way, Klein-Gordon equation has no source 
term: 



□ $ - m 2 $ = 

Equation (^) gives the scalar-field (bosons) energy-momentum tensor: 

t£ = ft,**0„$ + d^*d v $ - g^{d x <S>*d x <S> + ™ 2 $*$) 
From the Nother theorem the action derived from (^) gives rise to a conserved current: 

and the conserved charge: 

- 9 J° 



d 3 ry 



(4) 
(5) 
(6) 

(7) 



which can be identified with the number of bosons. 

Once the scalar field describes particles a quantum analysis is required. If we take Nb bosons in their lowest state, 
| Nb), the Einstein equations G M „ = (Nb \- T^ : \Nb) give exactly the same set of equations. This result is obtained 
with the scalar field $ changed into the operator: 



$(r, t) 



E 

k 



f 



ak0fc(r)e" 



(8) 



where creates quanta associated with the field <E> and au destroys them, obeying the commutation relation 
Sij . Considering the ground-state, one can introduce a semiclassical scalar field: 



'N B + l/2 



<M0 e ~ 



(9) 



where </>o and ljq are the eingenfunction of the corresponding lowest eingenenergy, which gives a nodeless eingenfunction 
4>o(r). Thence, forthcoming equations and numerical results are related to the semiclassical field $ c . 
The requirement of the functions to be orthonormal: 



allows one to determine the number of bosons, Nb, by using the classical field 



t> UJQ 



(10) 



(11) 
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B. Introducing Fermions 



Boson- fermion model star has been originally proposed by Henriques et al S, and are equilibrium compact cold 
configurations made from both bosons and fermions. 

As proposed by Chandrasekhar j(| for neutron star models, an equation of state can describe a perfect fluid of 
degenerate Fermi gas, and density of energy and pressure are given by: 

p = K(sxnht - t) 

(12) 

K ( . , . t 

p = — I sinh t — 8 sinh — — 'it 

where K — m*/327r 2 , and m n is the fermion mass (which is called neutron mass for neutron stars). The parameter t 
is related to the maximum momentum value, Qq = qo(j), in the Fermi distribution in a distance r from the centre: 

(13) 




The density of fermion particles in the distribution is: 



^) = ^ = ^§sinh 3 ( 7 ) (14) 



3^3^ Smh V 

Considering a compact spherically symmetric distribution with ([!]) and (|lj), the number of fermions is: 

iV F = J 47rr 2 nVAdr (15) 

Oppenheimer and Volkoff jl| have found equilibrium configurations by using Einstein equations with T± — Tf = 
T3 = p and Tq = —p, and with the hydrostatic equilibrium equation for the pressure: 

,' = -^f (16) 

where ' = d/dr, and B — B(r) is defined on (|l]). 

As proposed by Henriques et al § , the basic properties of the composed boson- fermion star can be studied with 
the use of an energy momentum tensor given by the sum of bosonic (B) and fermionic (F) counterparts: 



where: 



— T^ IV + (17) 



T lu = {P + P)u^,u v + g^p (19) 

where m m = (u°, u r , u$, u v ) = ^f- is the fermion fluid four vector. 

The set of differential equations can be obtained by using Einstein equations together with (|]) and (|l6|). After 
some calculations and by using variables redefinitions: 

x = mr , <y{x) = \f&nG<p ( r ) 1 w = — 

m 

(20) 
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we obtain: 



A' = zA 2 





/w 2 




<t' 2 " 


2p + 




m) <t 2 h 


A 



(21) 



5' = xAB 



2p 



w 



i U z 



.4 



(A-l) 



(22) 
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er' - A 










4); 

















(23) 



B' sinht - 2sinh(t/2) 
'1? cosh£-4cosh(i/2) + 3 



(24) 



Also, any function -F(r) is replaced with F{x/m). Hence, F'(x) = ^E^hul = _L dF ^ , where now ' = d/dx. 
Equations ()2l]) - J22] ) are obtained directly from Einstein equations using metric (0), and (||) is the Klein-Gordon 
equation (IJ) - Equation (^4|) arises from hydrostatic equilibrium equation (|l6|), rewritten in the form: 



f' 



1 g 

dpJdtW 



(P + P) 



(25) 



where: 



p= #(sinht-i) 



A' 



(sinhi- 8sinh(t/2) + 3i) 



(26) 



with: 



K = 



%imi 2 M 2 x 



(27) 



For convenience we will use K = 1/4-7T, which is the same value used by Oppenheimer and Volkoff for neutron stars. 
Note that (|27]) displays a constraint between boson mass m and fermion mass m n . Hence, if one choose m n as the 
neutron mass and K = l/4ir, we are setting m = 5, 11 x 10~ 17 MeV. 



III. SLOWLY ROTATING BOSON-FERMION STAR 

Slowly rotating relativistic stars has been studied for long in applications with neutron stars p0| . Meanwhile, until 
Kobayashi et al |ll|| rotation configurations for boson stars has never been studied, and until Silveira et al [fl2"| stable 
boson stars rotating configurations has never been found. Indeed, Kobayashi et al show that there is no stable config- 
uration with nonvanishing angular momentum for pure boson star studied perturbatively. Meanwhile, considering the 
quantum background for an axisymmetric scalar field one obtains stable configurations in the Newtonian framework. 

In this section we study relativistic equations for slowly rotating boson-fermion stars, in a similar framework 
presented by Kobayashi et al. 

A. Rotating Fermion Stars 

Using Hartle's prescription |l(J for slowly rotating neutron stars, one introduces the metric: 

ds 2 = -H 2 dr 2 + Q 2 dr 2 + r 2 K 2 \d9 2 + sin 2 9 {dip - Ldrf] (28) 
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where H, Q, K and L are functions of r and 9. When considering slow rotation we have Rfl <§C c, where R is the 
average radius and f2 is the angular velocity as seen by an observer at infinity. In this section we expand the metric 
( p8| ) only up to order fi 2 . This way, energy and pressure can be approximated using: 

s = P + o{n 2 ) 



(29) 



v = P + o(n 2 



When rotation is considered this metric introduces a dragging of inertial frames. Thence, L = dip/dr is the angular 
velocity acquired by one observer free falling the infinity to the point (r, 9). 

The metric ( p8|) must be invariant under rotation reversion, i.e., tp — > — ip and ft — ► —ft. Thus, functions H, Q and 
K are even in powers of Q; meanwhile L is odd. Since for slow rotation we study terms only up to order fl 2 , one can 
consider effects up to order f2 in L(r, 9), represented by C(r, 9): 

L(r,9) = C(r,9) + 0(n 3 ) (30) 

In cases when there is no rotation, Q = and the energy- momentum tensor for the perfect fluid is: 

T£ = (p + p)u»ti v +p5i; (31) 

But, when fl ^ 0: 

Tl l = {£ + V)u»u v + P5£ (32) 

In the case analysed by Hartle |lC|| , the first order term in (|30|), C(r, 6*), can be solved analytically, using the 
equation: 

R° 3 = 8irGT° (33) 
we are considering it M = (u T , u r , u e , u v ^j = (ti^O, 0,u 3 ) due to symmetries involved. Since dip = ildr, one can use 



u 
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Qu°. With the normalization condiction u^u^ = — 1, one obtains [<?oo + 2£!go3 + ^333] = — 1- Thus, the 

energy- momentum tensor component T® is: 

T 3 ° = (£ + T)u 3 u° = (£ + V) [(u°) 2 g 03 + u°u 3 g 33 ] = (£ + V)(u°) 2 (g 03 + n 933 ) (34) 

Expanding this term up to order fl 2 yields: 

^3° = (p + p)e~ u (n-C)r 2 sm 2 9 + 0(n 3 ) (35) 

Therefore, the field equation becomes: 



r 4 



,(X- u )/2 



r 4 e -(^)/2 drC + e - _ g e f sin 3 gdgC \ = iQ 7T G(p + p)e^- l, ^ 2 C (36) 



where C — £1 — C '. Defining the quantity: 

j(r) = e -(-+ A )/ 2 (37) 

one obtains: 



dj 



dr 2 

From Einstein equations: 



+ \')e-^l 2 (38) 



A' = %TtGre x p- -(e x - 1) 
r 
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(39) 



v' = 8irGre x p+ -(e A - 1) 
r 



one is able to rewrite (pq ) as: 



r dr 



Thus (pq) becomes, up to first order in fi: 



\d r [r*j(r)d r C] + -(d rJ )C + ^ , V/ l de (sm 3 8d e C) = 



Expanding C in vector spherical harmonics: 

oo 

and thus: 



z=i 



smf 



-d e Pi(cos6) 



-d r [r^j^drd] + d 



= 



(40) 



(41) 



(42) 



(43) 



When r — > oo, A and ^ vanish, and j —* 1. Thence, when r — > oo in Q43J), C; behaves like: 

C;(r) -> fcir - ' -2 + fc 2 r /_1 (r -> oo) (44) 

where fci and fc2 are constants. Since fields decrease very fast the solution when r — > oo is of the Kerr-Newman form: 

C; oc 

in such way that Ci = Vl — Q ~ Q. This is reproduced by (^i|) by choosing 2 = 1, and appropriate values for fci and 
fc 2 . Taking again correct values for the constants in such wise the solution to be regular, one can observe that terms 
with I > 1 decrease even faster and I = 1 mode can be considered as the dominant one. Hence, C(r, 9) ~ C%(r), and 
C becomes function of r alone. Now, (43) is: 



\d r [r 4 j(r)9 r d] + -{d r j)Ci = 



(45) 



Since the solution must be regular close to the origin, Ci(0) =const., and considering the exterior solution, one 
obtains: 



Ci{r) = n- 2 ± 



(46) 



where J is the total angular momentum of the fermion star. The quantity J/tt defines its moment of inertia. 



B. Introducing Bosons 

After reviewing Hartle's perturbative approach for rotating fermion stars (which are usually referred as neutron 
stars) it is convenient to introduce bosons to obtain the composite boson-fermion star model, which is the original 
part of this article. 

Following the perturbative approach used by Hartle, we expand the metric up to terms in Q 2 , and the Kerr metric 
with no charge can be written as: 

ds 2 = -B(r)dr 2 + A(r)dr 2 + r 2 d6 2 + r 2 sin 2 ^ 2 - 2C{r)r 2 sin 2 edrdip + 0(n 2 ) (47) 



G 



where C (r) = C\ (r) in the previous section. The total energy- momentum tensor is again: 



rp rpa I rph 

as in equation (frlj). Up to order f2 2 , <p= f2, and the nonvanishing components of T M „ are: 



J 



(48) 



(49) 



us 
~B 



il 

A 



(50) 



n 



B 



A 



(51) 



n0 = (P + P) r 2 sin 2 0[r! ^ C] 



(52) 



T 3 = (p + p)f! 



(53) 



Using redefinitions (|2(]) and after some straightforward calculations, one obtains up to order SI 2 , the set of differential 
equations: 



A' = xA 2 



2p+ { — ■ }^ 



T ' 2 



A 



(54) 



B' = xAB 



2p 



w 
~~B 



I ff 31 



_/ 2 "I 



.4 



(A-l) 



(55) 



1 

2 IB 



(56) 



t' = -2 



B' sinht - 2sinh(i/2) 
~B coshi-4cosh(i/2) + 3 



(57) 



C" 



"4 








.r 
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C -AA{p + p) (O-C(r)) 



(58) 



where now ' = d/dx. 

Equations (|54|) to (^7|) are the same as those obtained for boson- fermion stars. Equation (|5^) gives the function 
C(x), corresponding to rotation of the boson- fermion star. Note that C(x) depends not only on x, but the parameter 
t and metric fields also contributes to the rotation. Thus, step-by-step results that numerically arise from ( f34|) - (|57| ) 
can modify the rotational term C(x), via changes on the initial values, Co ancl to, which represents bosons and fermions 
contributions respectively. On the other hand, C{x) does not appear on equations (54) - (p7|). Hence, changes on the 



initial value of C(x) does not produce changes on cr(x), t{x), A{x) or -B(cc), up to order fr 
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IV. RESULTS 



After obtaining perturbatively, up to order fi 2 , the equations that govern the rotation of a boson- fermion star, the 
results from section 3.2 are obtained numerically. Before that, equations (|5^) and ( [58] ) ought be rewritten to form a 
first order set of differential coupled equations. Hence, equation (^) splits into: 

s = a' (59) 



S ' = (xAo 2 ±ti) s - Aa - l) (60) 

and taking: 

C(x) = CI - C(x) (61) 

equation (|||) becomes: 

C' = -U (62) 



U - AA (p + p) C (63) 

The boundary condictions are Ao = 1, sq = 0, Uq = U(Q) = 0, considering Bo, ao, to and C(0) = Co as arbitrary 
initial values to be ascribed. And, at the infinity condictions are A(oo) — > 1, B(oo) — > 1, ct(oo) — * and s(oo) — * 0. 
As seem before, star properties as mass, radius and total number of particles does not differ from boson-fermion 
stars with no rotation up to fl 2 . For example, a configuration with ao = 0.2, Co = 1 and to = 4.0 we have obtained 
M ~ QAlMpJm, (raN-B +m n ATp) ~ OAQMpJm e mR ~ 0.94M|j/m. Thus, if one finds an equilibrium configuration, 
changes on rotation parameter C does not modify those results. Meanwhile, if one fix the initial value Co, changes on 
values of ao an d ^o provide modifications on the C values evolution, i.e., the referential frames dragging is modified. 

Evolution of the scalar fields, its derivative and metric coefficients are shown in figure 1, and rotation parameters 
are shown in figure 2. Those figures are obtained with ao — 0.20, to = 4.0 and Co = 1.0, for which Bo — 0.1422 
and w = 0.79108. In this case, E-q < 0, showing that configuration is probably stable. We can also observe that C 
stabilizes for x — > oo. But this is expected and from equation (|46| ) one can see that this maximum value is Q, up 
to the approximation order we assume. One can obtain this by using the boundary condictions at infinity, where 
asymptotic flatness is expected for (p8|), and lim r ^oo A(r) = 1 = limj-^oo B(r) and lim r ^oo C(r) = 0. The last one is 
equivalent to: 

lim C(x) = (64) 

x — >oo 

Since C(x) = f2 — C(x), equation ( |64[ ) is possible if and only if: 

lim C{x) = Q (65) 

x — >oo 

In the case shown in the figures 1 and 2, angular momentum as seem by an observer far from the object is 
ft = —0.1651. Some other results are shown in the table I below for a configuration with to = 4.0 and ao = 0.2. 



U' 



xA 



P+ ^ a ~A 



4 

x 



Co 


3.0 


2.0 


1.0 


-1.0 


-2.0 


Q 


-0.4954 


-0.3302 


-0.1651 


0.1651 


0.3302 



TABLE I. Values obtained for il given an initial value of the rotational parameter Co, with ao = 0.2 e to = 4.0. 
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All configurations obtained present £b < 0, and the same values for mass and radius. Note that negative values 
just provide an inversion of rotation, as expected. 
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FIG. 1. Scalar field, its derivative and metric coeffcients for a composite boson-fermion star, with <to 
Co = 1.0. Values obtained: B = 0.1422, w = 0.79108 and E B ~ -0.05. 



0.20, t = 4.0 and 




FIG. 2. Rotation parameter evolution for the case shown in figure 1. In this case lim^^cx) C(x) — £1 = —0.1651. 



V. CONCLUSIONS 



Given a star composed of both bosons and fermions, without interaction, we have show that it is possible to obtain 
stable configurations with slow rotation, by using the perturbative rclativistic method that usually describes neutron 
stars. We have developed sample calculations considering fields up to the second order in angular velocity. With that, 
angular velocity, mass, number of particles and binding energy can be measured by an observer at the infinity. 

But, if it is known ]ll| that the perturbative prescription allows no rotation for boson stars, why have we obtained 
rotation with the same prescription when bosons are mixed to fermions inside the star? The answer is in equation 
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(fig). Up to the order considered here, rotation is not directly influenced by the scalar field; instead, energy density 
and pressure of the fermi gas are explicit in the equation that describes the rotation of the configuration. 

Since we have performed I = case, which is the predominant one, these results may contribute to studies on 
deformations of the structure, e.g. I — 2, and so on, and studies on the emission of gravitational waves from those 
kind of objects. This study also opens some questions about spinning boson fermion stars with the use of axisymmetric 
scalar field. 
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